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ABSTRACT 

We extract J and K s magnitudes from the 2MASS Point Source Catalog for approx- 
imately 6xl0 6 stars with 8 < K s < 13 in order to build an Ak extinction map 
within 10° of the Galactic centre. The extinction was determined by fitting the up- 
per giant branch of (K s , J-K s ) colour-magnitude diagrams to a dereddened upper 
giant branch mean locus built from previously studied Bulge fields. The extinction 
values vary from =0.05 in the edges of the map up to Ak =3.2 close to the Galac- 
tic centre. The 2MASS extinction map was compared to that recently derived from 
DENIS data. Both maps agree very well up to ^4^- =1.0. Above this limit, the compar- 
ison is affected by increased internal errors in both extinction determination methods. 
The 2MASS extinction values were also compared to those obtained from dust emis- 
sion in the far infrared using DIRBE/IRAS. Several systematic effects likely to bias 
this comparison were addressed, including the presence of dust on the background of 
the bulk of 2MASS stars used in the extinction determination. For the region with 
3° < |&| < 5°, where the dust contribution on the far side of the Galaxy is ~ 5 %, 
the two extinction determinations correlate well, but the dust emission Ak values are 
systematically higher than those from 2MASS. A calibration correction factor of 76% 
for the DIRBE/IRAS dust emission extinction is needed to eliminate this systematic 
effect. Similar comparisons were also carried out for the 1° < |&| < 3° and \b\ < 0.5° 
strips, revealing an increasing complexity in the relation between the two extinction 
values. Discrepancies are explained in terms of the calibration factor, increasing back- 
ground dust contribution, temperature effects influencing the dust emission extinction 
and limitations in the 2MASS extinction determination in very high extinction regions 
(| &| < 0.5°). An asymmetry relative to the Galactic plane is observed in the dust maps, 
roughly in the sense that Ak values are 60% smaller in the south than in the north 
for 1° < |&| < 5°. This asymmetry is due to the presence of foreground dust clouds 
mostly in the northern region of the Bulge. 
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1 INTRODUCTION 

The high extinction and its patchy distribution in the Galac- 
tic centre region have been a constant problem to the study 
of the properties of the Bulge stellar population. Several 
efforts were carried out to investigate the extinction distri- 
bution close to the Galactic centre. For example, Catchpole 
et al. (1990) studied the distribution of stars in the cen- 
tral 1° x 2° of the Galaxy by means of J, H and K colour- 
magnitude diagrams (CMDs). They derived visual absorp- 
tions in the range 7 < Av < 30. Frogel et al. (1999, hereafter 
FTK99) obtained extinction values varying from Ay=lA\ 
up to yl v =19.20 for 11 Bulge fields with |6| < 4°, close to 
the Galactic centre. Stanek (1996) derived a mean extinc- 
tion of < A v > = 1-54 for Baade's Window (Baade 1963). 
As a consequence, the investigations about the Bulge stellar 



content based on optical data were restricted for a long time 
to lower-extinction regions such as Baade's Window. 

The advent of near infrared surveys such as the Two 
Micron All Sky Survey (2MASS, Skrutskie et al. 1997) and 
the Deep NIR Southern Sky Survey (DENIS, Epchtein et al. 
1997) has provided fundamental tools to study the stellar 
population (Unavane et al. 1998) and extinction (Schultheis 
et al. 1999) in the inner Bulge. Recently, Dutra et al. (2002, 
hereafter Paper I) confirmed the existence of two new Bulge 
windows, W0. 2-2.1 and W359.4-3.1, closer to the Galactic 
centre than Baade's Window; they used the JK S photometry 
from the 2MASS survey archive to map the extinction dis- 
tribution within 1° towards these windows. The reddening 
distribution in the Bulge area is also fundamental to under- 
stand the spatial distribution of globular clusters and their 
relation to the Bulge field stellar population itself. Barbuy et 
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Figure 1. Extinction-corrected (filled circles) and observed (open 
circles) CMDs for a cell at t = 1.93°, b = 0.93°. The straight line 
represents the reference upper giant branch (Eq. 1). 

al. (1998) discussed the globular clusters projected within 5° 
of the nucleus and which appear to be related to the Bulge, 
while Barbuy et al. (1999) discussed the properties of those 
found in the area covered by the present study. 

In this work we use the 2MASS J (1.25^im) and K s 
(2.17/um) photometric data to build (K s , J-K 3 ) CMDs of 
Bulge fields in order to map out the interstellar extinction 
in the central 10° of the Galaxy. The extinction was derived 
by means of the upper giant branch fitting method, similar 
to that described in Paper I. In Sect. 2 we revise this method 
and describe in detail the way we build the extinction map. 
In Sect. 3 we compare the present results with those from 
DENIS. In Sect. 4 we analyze the 2MASS and DIRBE/IRAS 
extinction maps, making use of simple models for the opti- 
cal depth in different directions in the Galaxy, and discuss 
the asymmetries between observations in the northern and 
southern Galactic hemispheres. Finally, the concluding re- 
marks are given in Sect. 5. 



2 BUILDING THE EXTINCTION MAP WITH 
2MASS DATA 

In Paper I we have built an extinction map in regions of 
low-extinction by means of upper giant branch fitting. We 
adopted an upper giant branch template from a composite 
CMD using 7 Bulge fields from FTK99. Using this template 
we managed to reproduce the mean extinction in Baade's 
and Sgr I windows. Metallicity effects were not considered 
in Paper I, since Ramirez et al. (2000), from a spectroscopic 
study of central M giant stars, pointed out that there is no 
evidence for metallicity gradient within the inner Bulge (R 
< 560 pc). Besides, FTK99 concluded that the amplitude of 
metallicity variations in the inner Bulge implies very small 
giant branch slope changes. Therefore, we assume that the 
metallicity variations in the inner Bulge do not affect sig- 
nificantly the extinction estimates and use as reference the 
same upper giant branch adopted in Paper I: 

(K s ) = -7.81(J - K B ) + 17.83 (1) 



Figure 2. Analysis of limiting magnitude for fields with \b\ < 
1.5°: CMDs for stars with (a) 8 < K s < 11.5 and (b) 9 < K s < 
13 within a cell at I = 0.87°, b = 0.27°; (c) and (d) same K s 
intervals as in (a) and (b), respectively, for stars within a cell at 
I = 1.93°, b = 0.93°. The straight line represents the reddened 
reference upper giant branch. The derived extinction values in 
panels (a) and (b) are A^= 1.88 and 1.83, respectively. For (c) 
and (d) A K = 0.74. 



The equation above describes appropriately the upper 
giant branch locus for stars with 8 < K o <12.5. 
We carried out JK S photometric extractions of stars in the 
2MASS Point Source Catalog available in the Web Inter- 
face http://irsa.ipac.caltech.edu/applications/Gator/ for 61 
fields with radius r^l° each. For most fields, we extracted 
stars in the range 8 < K s < 11.5, there being typically 60,000 
such stars in each. For the fields with \b\ < 1.5°, we ex- 
pected very high extinction (Schultheis et al. 1999); thus 
we extracted fainter stars, down to K s =13. This extra 1.5 
mag allowed us to determine the best magnitude range for 
upper giant branch fitting, taking into account factors such 
as the increasing photometric errors and contamination by 
disk stars with fainter magnitude limits, and the decreas- 
ing observed extent of the upper giant branch in heavily 
reddened inner Bulge fields. This issue is discussed in Sect. 
2.1. In these low-latitude regions the deeper 2MASS extrac- 
tions led to a considerably larger number of stars per field, 
ps 200,000. Considering all fields in the area, we extracted J 
and K a magnitudes for approximately 6xl0 6 stars. 

In order to map out the extinction in the 10° Galactic 
central region, we defined 961 cells with 4' x 4' in size in 
each one of 61 extracted fields. The extinction was deter- 
mined assuming that the upper giant branch defined by the 
stars in a cell has the same slope as that adopted as reference 
(Equation 1). We then calculated from the (K, J-K s ) values 
of each star the shift along the reddening vector necessary to 
make it fall onto the reference upper giant branch. We used 
the relations A K , = 0.670£( J - K„) and -^f - = 0.95 (Paper 
I) to derive the Ak values. The Ak adopted for each cell 
was taken to be the median of the distribution of Ak values 
for the individual stars in it; an iterative 2-a clipping was 
applied to the distribution of Ak values in order to elimi- 
nate contamination from foreground stars. Fig. 1 shows the 
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Figure 3. Ak,2MASS extinction map within 10° X 10° of the Galactic centre. The contours represent the following extinction values: 
A K = 0.2, 0.5, 1.0, 1.5 and 2.0. 



method applied to a cell located at I = 1.93°, b — 0.93°; 
the stars in the observed CMD (open circles) have a median 
extinction of Ak = 0.74 with respect to the reference upper 
giant branch. The extinction-corrected CMD based on this 
value is also shown (filled circles) in the figure. Note that 
the extinction-corrected CMD is well fitted by the reference 
upper giant branch (straight line), indicating that the de- 
termined median extinction is representative of most stars 
within the cell's area. 



2.1 The upper giant branch fitting in heavily 
reddened fields 

We now try to assess the existence of systematic effects on 
the extinction determination just outlined. Pushing the lim- 
its used in our fitting method towards fainter magnitudes 
is desirable only if it increases the magnitude range where 
the upper giant branch is clearly defined and can be fitted. 
Several factors may prevent fainter magnitudes from being 
useful in this way. The first is obviously the increasingly 
large photometric errors. The 2MASS JK S nominal errors 
are smaller than 0.04 for J ~ 15 and K s ~ 11.5. Close 
to the centre, however, the density of stars increases very 
rapidly, yielding less precise photometric measurements due 



to source crowding. Another issue is the existence of selec- 
tion effects: we fit a straight line to the upper giant branch 
defined in a (K s , J — K s ) CMD. The K s limit to the data 
will stand out clearly as a lower limit in the populated CMD 
region, but any J band cut-off will be a diagonal line on the 
CMD. It will bias our best-fit line if the giant branch is 
crossed by this diagonal line within the magnitude range 
used in the fit. Finally, at fainter magnitudes the Bulge lu- 
minosity function rises rapidly, increasing the contamination 
by Bulge stars with less extinction than the average, which 
will hamper our attempt to fit the bulk of the Bulge stellar 
population. 

All the 3 factors just mentioned are in fact present in our 
data to some extent. Fig. 2 shows the CMDs for two sample 
cells: one at £ = 0.87°, b — 0.27° (panels a and b), the other 
at I = 1.93°, b = 0.93° (panels c and d). The panels on 
the left (right) show the data in the range 8 < K s < 11.5 (9 
< K s < 13). The reference upper giant branch is shown in all 
panels, reddened by the best-fit (Ak,E(J — K)) value found 
in each case. Clearly, the extinction is underestimated for the 
lower latitude cell when stars in the 9 < K„ < 13 interval 
are used. The reason for that is two-fold: there is a large 
number of blue stars with K s ~ 11 and there is a paucity 
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Figure 4. Histogram of Ak extinction values of the 10° extinc- 
tion map. 

of stars in the lower right corner of the CMD in panel 2b. 
Both tend to shift the best fit line towards the blue. 

The scarcity of faint red stars is due to the 2MASS J 
band limit: there are essentially no stars beyond the diag- 
onal line J = 16.5, which is thus the 2MASS Point Source 
Catalog limiting J magnitude. This empirical limit is close 
to that proposed by the 2MASS collaboration (Skrutskie et 
al. 1997, see also the 2MASS web site). The faint blue stars 
are probably due to contamination from Bulge stars with 
less extinction and to the increased photometric errors in 
these dense fields. 

Note that the fit within the range 8 < K s < 11.5, 
despite the inclusion of some residual stars bluewards of 
the giant branch is not biased (panel 2a). For the cell at 
£ = 1.93°, b — 0.93° the situation is more reassuring, the 
derived Ak being insensitive to the magnitude range used. 
We thus conclude that a more efficient use of the CMDs in 
the crowded and high extinction areas close to the Galactic 
plane is made by restricting the fit to K S =11.Q. This limit is 
also conveniently close to where completeness effects should 
start to be significant in the 2MASS. For DENIS, Unavane 
et al. (1998) estimated the 80% completeness level to be 
at K ~ 10; J ~ 13 and the 2MASS data are about 1 mag 
fainter. 

We conclude that in areas where extinction is around 
Ak < f .5 and photometric errors are not substantially in- 
creased by crowding our fitting method does not suffer from 
any biases due to the J band cut-off. For Ak ^ 2.5, the J 
band detection limit effects become dominant and the range 
available for the fit is substantially shortened, rendering the 
extinction determination unreliable. 

2.2 The Ak extinction map 

Fig. 3 shows the Ak contour map for the central 10° x 10° 
of the Galaxy obtained by applying the method described in 
the previous section to stars down to K S = 11.G. The Ak con- 
tours show that the regions with Ak > 1.5 are concentrated 
in the area with \b\ < 1.0°. The isocontours with high val- 
ues are slightly asymmetric with respect to the mid-plane 
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Figure 5. Upper panel: histogram of internal error di values. 
Lower panel: Variation of Ui with Ak extinction values. 



and more extended in the southern hemisphere. They ap- 
pear to be caused by the displaced location of the Sun from 
the Plane (Sect. 4). The structure with extinction values be- 
tween Ak—0. 5-1.0 and angular dimension of 3° x 2° around 
I = 1.5° and b = 4.0° is a component of the Pipe Nebula, a 
dark nebula recently studied in CO by Onishi et al. (1999). 
This nebula appears to be located at 160 pc from the Sun, as 
a southern extension of the Ophiuchus dark cloud complex, 
and is located at the edge of the ScoOB2 association (On- 
ishi et al. 1999). The detection of the Pipe Nebula in Fig. 3 
suggests that nearby clouds cannot be neglected in attempts 
to interpret central extinction maps. From the 32,761 cells 
that cover a 12° x 12° central area of the Galaxy, we have 
2MASS data for 80 % of them. The quadrant 0° < £ < 5°, 
—5° < b < 0°, which comprises Baade's Window, is so far 
only partially released by 2MASS. 

Fig. 4 shows a histogram of A K values for the cells in the 
10° x 10° extinction map. The mean extinction in the entire 
map is < Ak > = 0.29 with a standard deviation a = 0.12 
from the mean. 63 % of the cells fall within 2-a of this mean 
value, and 80% of them have Ak < 1-0. The upper panel of 
Fig. 5 shows the histogram of internal errors in extinction 
determination. The mean internal error is < o» >= 0.08, 
with a standard deviation of 0.02 around this mean. 70% of 
the cells have internal errors within 2 standard deviations 
from the mean value (0.04 < Oi < 0.12). The lower panel 
shows the dependence of internal errors with Ak; we note 
that for Ak > 1.5 the internal errors increase significantly. 

Since the present extinction map can be useful for a 
wide variety of Galactic and extragalactic studies in such 
central directions, it will be provided in electronic form in 
the CDS, by columns: (1) and (2) galactic longitude and 
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Figure 6. Comparison between Ak ,2MASS and Ak , DENIS extinc- 
tion values for the region \£\ < 5° and \b\ < 1.5°. The straight 
line represents the identity function. 



latitude of the cell centre, (3) the K band extinction Ak 
and (4) the uncertainty in the Ak determination Uj. 



3 A 2MASS VS. DENIS EXTINCTION 
COMPARISON 

Schultheis et al. (1999) provided an extinction map for the 
inner Galactic Bulge (covering \£\ < 8° and \b\ < 1.5°) 
obtained from the J and K s DENIS CMDs, together with 
isochrones from Bertelli et al. (1994), hereafter called DENIS 
extinction map. The differences between the present extinc- 
tion determination applied to 2MASS photometry and the 
procedure adopted by Schultheis et al. (1999) are the fol- 
lowing: (i) the latter use an isochrone from Bertelli et al. 
(1994), with metallicity Z—0.02, age of 10 Gyr and distance 
d=8 Kpc, to represent the dereddened (K,J-K) CMD of the 
Bulge fields; (ii) Schultheis et al. adopt a transformation 
from K to K s with an estimated error of 0.04 mag, in or- 
der to be able to use isochrone fitting; and (iii) the DENIS 
infrared photometry is limited to 7v" s = 11.0, with detection 
limits at J = 16.0 and K s = 13.0 (Schultheis et al. 1999). 

Unavane et al. (1998) estimated a DENIS completeness 
limiting magnitude of K a = lQ.Q in the inner Bulge. The res- 
olution of the DENIS extinction map is also 4'. Fig. 6 shows 
the comparison of the extinction values derived from the 
2MASS photometry, Ak ,&mass, with those derived from the 
DENIS photometry, Ak,denis, for the area in common be- 
tween the two maps (\£\ < 5° and \b\ < 1.5°). The ex- 
tinction values derived from the 2MASS and DENIS pho- 
tometric data present an excellent agreement, especially up 
to Ak=1'0. Beyond this limit Ak.denis values are higher 
than Ak,2mass ones, but if we consider that the uncertain- 
ties in the extinction determination and photometric errors 
increase in these zones, the agreement is still significant. The 
departure from identity line in Fig. 6 could be partially due 
to differences in the filter profiles and zero-point calibrations 
adopted by the two surveys. 
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Figure 7. Increase in optical depth as a function of distance from 
the Sun as predicted by a simple model for the dust distribution 
in the Galaxy. The upper panel shows directions along the disk 
plane, for two values of longitude I, as indicated. The mid (lower) 
panel shows directions with 6 = 2° (4°). 



4 A 2MASS VS. FIR EXTINCTION 
COMPARISON 

Schlegel et al. (1998) (hereafter SFD98) presented an all-sky 
reddening map based on the 100 /im dust emission, model- 
ing the emission by dust grains with blackbody radiation 
at a temperature T=48.2 K. Temperature corrections were 
applied by means of the 100 and 240 /im DIRBE maps, 
up to T = 21K and in low angular resolution (1°). Their 
E(B-V)fir values correspond to the integrated dust col- 
umn contribution throughout the Galaxy. The transforma- 
tion from E(B-V)fir to Ak,fir assumed Ak = 0.112^4^ 
and R v = A V /E(B — V)= 3.1 (Cardelli et al. 1989). The 
2MASS extinction values, on the other hand, are based on 
the position of Bulge red giant stars in the CMD and, as 
such, are limited in depth by opacity effects in the employed 
bands, especially in the J band. 

In this section we compare the Ak extinction values de- 
rived from dust emission, Ak,fir, to those from the 2MASS 
photometry, Ak.2mass- But before we can properly com- 
pare these two extinction maps, we must first investigate the 
possible existence of systematic effects on this comparison, 
such as variations in dust temperature and line of sight dis- 
tribution or the presence of high density intervening clouds. 
Such effects should arise due to the widely different obser- 
vational signatures from the interstellar medium on which 
the two extinction values are based. 
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Figure 8. Relative loss of sensitivity in the 2MASS extinction de- 
termination method due to high density intervening dust clouds. 
The optical depths are assumed to be averages over the entire 
cell of 6.1' in side. The cloud column density profile is assumed 
to be exponential, with scale length 6>o and central column den- 
sity To - Tpm is averaged over the entire cell but the 2MASS value 
is assumed to be restricted only to points with t < 2. The ratio 
t 2MASs/ t fir. is plotted against to for different values of 6o, as 
indicated. 



4.1 Foreground dust vs. background dust 

In order to help interpreting the relation between Ak,2mass 
and Ak,fir, we consider a simple model of dust distributed 
exponentially along and perpendicular to the Galactic plane. 
The optical depth out to some distance r from the Sun in 
the direction given by Galactic coordinates (£,b) is then: 



r (r, I, b) oc 



Jo 



-R(r,l,b)/R d -\Z(r,b)\/Z d 



(2) 



where 

Z — rsinb and R 2 — Rf, + r 2 cos 2 b — 2 Ro r cosb cost are 
cylindrical coordinates centred on the Galaxy. Rd and Zd are 
the dust horizontal and vertical exponential scales, whose 
values we assume to be 2.5 kpc (Robin et al. 1996 and 
Drimmel & Spergel 2001) and 110 pc (Mendez & van Al- 
tena 1998), respectively. We also adopt Ro = 8.5 kpc for 
the Sun's distance to the centre of the Galaxy. We numer- 
ically performed the integral above for the same directions 
for which we have 2MASS and DIRBE/IRAS data. For each 
direction, we assume that our model Ak,2MASS and Ak,fir 
values are proportional to t(Ro) (hereafter rg.s) and t(oo) 
(hereafter Too), respectively. To model the foreground and 
background distribution with respect to the Galactic Cen- 
tre, we also assume that the solar position is displaced by 
Zaun = 15 pc above the Galactic plane (e.g. Cohen 1995); in 
the double exponential dust distribution model, the fraction 



of forcground-to-total dust distribution does not depend on 
whether the Sun is displaced above or below the Galactic 
Plane. A full data vs. models comparison is currently under 
way. For now we restrict our discussion to the basic model 
features and their differences relative to the data. 

Fig. 7 shows t as a function of distance from us for 
several directions (I, b) as predicted by our simple model. 
The figure clearly shows the dependence on Galactic lati- 
tude of the expected contribution of dust beyond Ro = 8.5 
kpc. The fraction of the total optical depth caused by dust 
beyond the Galactic centre varies from 50 % at the Galactic 
Plane to only a few percent for b — 4°. Fig. 7 also shows 
that the total optical depth itself decreased by nearly two 
orders of magnitude between these two Galactic latitudes. 
The dependence on £, on the other hand, is quite small. 



4.2 Resolution elements and zones of avoidance 

As we determined Ak,2mass in cells of 4' on a side, the 
resolution of our maps is larger than the 6.1' resolution of 
the SFD98 extinction maps. In order to place both maps 
on a similar angular resolution, we convolve the Ak,2mass 
extinction map with a <j = 4.5' Gaussian. 

However, there are biases that do not depend on the 
resolution scale of the maps. One important issue is that the 
Ak,2mass values require a minimum number of stars along 
the upper giant branch to be determined. Regions with Ak 
~ 2 may be obscured enough that only the brightest stars 
will fall in the range K s < 11 used to determine Ak,2mass 
(see Sect. 2.1). If the region covers the entire cell, then no 
extinction value will result. But more common will be the 
cases when this zone of avoidance partially covers a cell. 
In that case, the cell may still have stars enough for the 
giant branch fitting method to be applied but the resulting 
Ak,2mass will be an underestimate of the true one. Note 
that the Ak,fir values will not suffer from such bias. We 
may model the latter simply as an average over the entire 
cell: 



TFIR = 



JcM T ^ )dU} 



where the integral is over the cell's solid angle and r(ui) is 
the varying K band optical depth within the cell boundaries. 
For the 2MASS, assuming that only regions with r < 2 will 
contribute with CMD stars, we would have: 



T2MASS = 



J <2 r(u)du 

J zell 



where the 2MASS integral limit on the numerator corre- 
sponds to the cell regions where r < 2. We may then es- 
timate the relative bias in the latter quantity by determin- 
ing the ratio T2mass/tfir for some model for t(u). An 
important case of interest, especially for the lowest b re- 
gions (|b| < 0.5° (see Fig. 3), will be that of a direction 
cutting through a dense dust cloud. We may consider that 
the cloud will dominate the dust column density in that di- 
rection and therefore ignore the contribution of foreground 
and background material. If we also assume the cloud to 
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Figure 9. Ak } 2MASs/Ak,fir ratio map within 10° X 10° of Figure 10. Apc^MASS (solid line) and Ax, fir (dashed line) 

the Galactic centre. Large (small) ratio values correspond to dark profiles in galactic longitude centred in the galactic latitudes: (a) 

(light) areas as indicated in the greyscale bar shown on the right. b = 0.48°, (b) b = 1.48°, (c) b = -0.48° and (d) b = —1.48°. 



have a central optical depth To, circular symmetry, and an 
exponentially falling dust column density, we will have: 

r(w) = t(6) = To e- e/e ° 

where 9 is the angle between the direction to the cloud centre 
and the direction considered and #o is the cloud exponential 
scale. 

Fig. 8 shows T2mass /tfir as a function of r for sev- 
eral choices of 9q. For 9o = 2' and To > 5, the entire cell 
would have r > 2 and therefore tk,2mass would become 
undetermined. We infer from the figure that the bias on 
the 2MASS data is relatively insensitive to the cloud pro- 
file shape but it is clearly dependent on profile normaliza- 
tion. For dust clouds with to < 5 (roughly corresponding 
to a central Av < 55), the systematic effect on Ak,2mass 
will be ~ 20%. Note that these ratios were computed as- 
suming the centre of the dust cloud to coincide with the 
cell's centre. For the more common off-center positions, 
A k ,2mass/A~k,fir ^ t k ,im ass I Tfir will be larger than 
shown in Fig. 8. Also, in the higher extinction regions there 
may possibly be more than one cloud per cell, thus increas- 
ing this systematic effect. 

4.3 Ak,2mass / Ak.fir ratio map 

Fig. 9 shows a gray scale map of the Ak,2mass/Ak,fir 
ratio, which provides information about the relative contri- 
bution of dust on the background of the Galactic centre, as 
background dust affects only Ak,fir- Since this is a ratio 
map between 2 observed extinction distributions, the result 



is not dependent on the calibration details. We note a global 
effect in the sense that Ak,2mass < Ak,fir, which will be 
discussed later. One clearly sees from the ratio map the ex- 
istence of a flat region of very low Ak,2mass/Ak,fir ra- 
tio. In this region closest to the Plane, the FIR extinction 
values reach much larger optical depths than the 2MASS 
data. Besides, the FIR extinction values are severely over- 
estimated by FIR emission of dust heated at temperatures 
above those assumed in the correction based on the 100 and 
240 /im DIRBE maps. There are also two symmetric dark 
regions close to the Galactic Centre, where Ak,2mass is 
comparable to Ak,fir- These structures suggest a conical 
region depleted of dust on either side of and centred at the 
nucleus. This interpretation is based on the more limited op- 
tical depth of the 2MASS extinction values, and on the fact 
that the features on the foreground of the Galactic centre 
are present both in the 2MASS and FIR maps and should 
therefore be canceled out. This conical structure is similar 
to those observed in emission line maps from some AGNs 
(e.g. Storchi-Bergmann & Bonatto 1991). It might be related 
to the central low luminosity AGN in our Galaxy, with an 
ionizing cone orientation nearly coincident with that of the 
Galaxy rotation axis. However, it may more readily be ex- 
plained by a more extended central dust disc hole producing 
a shadow effect on the ionizing radiation. It would deplete 
dust grains within a cone, and the UV radiation sources are 
more probably related to continuous star formation in the 
central parts (e.g. Arches and Quintuplet clusters presently) 
and/or Post-AGB stars from the central Bulge (Binette et 
al. 1994). 

Another way of assessing the effect of dust depletion in 
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Figure 11. Comparison between Ak,2MASS an d Ak,fir. ex- 
tinction values for cells with 5° > b > 3° (upper panels) 
and —3° > 6 > —5° (lower panels). Panels (a) and (c) show 
the Ak 2M ASS I Ak FIR histogram, whereas panels (b) and (d) 
show the comparison over the entire range in Ak- The dashed 
lines in the latter panels represent the identity function, whereas 
the solid lines have angular coefficients equal to the median 
Ak,2MASs/Ak,fir value within each strip. 



Figure 12. Comparison between Ak : 2MASS an d Ak,fir. ex- 
tinction values for cells with 3° > b > 1° (upper panels) and 
— 1° > b > —3° (lower panels). Panels (a) and (c) show the 
Ak 2MASs/Ak fir histogram, whereas panels (b) and (d) show 
the comparison up to Ak = 2.0. The dashed lines in the latter 
panels represent the identity function, whereas the solid lines have 
angular coefficients equal to the main peak in the corresponding 
A K ,2MASS/A K ,FIR histograms. 



the cone region is to show cuts along galactic longitude at 
different latitudes. This is done in Fig. 10, where we show 
profiles of Ak,2mass (solid line) and Ak ,fir (dashed line) 
along galactic longitude in positions towards the conical 
structure. We clearly see the effect of central dust depletion 
in panels (10b) and (lOd), which is not expected at these 
galactic latitudes. 

One could naively argue that Ak,2mass < Ak,fir be- 
cause of dust emission beyond the Galactic centre, which 
contributes only to Ak.fir- Since the relative contribution 
of background dust should be the major reason for any dis- 
crepancy between 2MASS and SFD98 extinction values, we 
will analyze them in separate for three regions, each one 
divided into a northern and a southern strip: (1) \l\ < 5° 
and 5° > |6| > 3°, (2) \£\ < 5° and 3° > |6| > 1°, and 
(3) |£| < 5° and |6| < 0.5°. This division allows assessment 
of the correlation between Ak,fir and Ak,bmass in regions 
with varying contributions by the background dust. We hope 
therefore to disentangle the effect on Ak,fir by dust on the 
far side of the Galaxy from other factors that may lead to 
discrepancies in the Ak.fir. vs. Ak,2mass relation. 

4.3.1 Strips 5° > |6| > 3° 

Fig. 11 shows the Ak ,fir vs. Ak,2mass relation in the region 
5° > \b\ > 3°. The upper panels correspond to cells in the 
northern Galactic strip (|£| < 5°, 5° > b > 3°), whereas the 



lower panels correspond to the southern (\£\ < 5°, —3° > 
b > —5°) strip. The Ak,2mass / Ak,fir ratio histograms in 
panels (11a) and (11c) show that the relative extinction val- 
ues have peaks at 0.75 (75 %) and 0.72 (72 %) in the north- 
ern and southern strips, respectively. These well-defined 
peaks imply a good linear correlation between the two ex- 
tinction values, as confirmed in panels (lib) and (lid). For 
the northern and southern strips we have the correlations 
Ak,fir = 1-33Ak,2mass and Ak,fir = 1.39.Ak, bmass, re- 
spectively, where the angular coefficients were derived from 
the median Ak,2mass/Ak,fir ratio values. In this region, 
the model predictions indicate that, for both northern and 
southern strips, the ratio of extinction on foreground of the 
Galactic to total extinction should be 96%, independent of 
the Sun's displacement from the disc mid-plane. From the 
median Ak,2mass/Ak,fir values we infer that a typical 
value for this ratio is » 73%. Therefore this discrepancy 
between model and observed extinction ratio suggests that 
the contribution from background dust cannot explain the 
observed A k ,2mass / A k ,fir ratio. 

Arce & Goodman (1999) also found a linear relation 
between the dust emission extinction and that derived from 
the stellar content in the Taurus Dark Cloud, with a slope 
that varies in the range 1.3-1.5. They attributed this dif- 
ference to the dust column density vs. reddening calibra- 
tion from Schlegel et al. (1998), which may yield an over- 
estimated Ay >0.5. Note that the slopes in our two strips 
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Figure 13. Comparison between Ak,2MASS an d Ak ,FIR extinc- 
tion values for cells with 0° < b < 0.5° (upper panels) and 
0° > b > —0.5° (lower panels). Panels (a) and (c) show the 
Ak,2MASs/Ak,fir histogram, whereas panels (b) and (d) show 
the comparison over the entire range in Ak- 



(with 2 > Ay > 10) are similar to those of Arce & Goodman 
(1999), which corroborates the idea of a Ak.fir calibration 
effect. 

We point out that assuming a significantly lower value 
for Ry would compress the Ak,fir scale, decreasing the 
differences between Ak,fir and Ak,2mass in Fig- H- How- 
ever, observational constraints on Rv do not support this 
possibility. Indeed, Gould et al. (2001) estimated Rvi = 
A V /E(V -I) = 2.4 (R v « 3.0), from VIK colours of 146 
Baade's Window G and K giants. Stanek (1996) obtained 
Rvi = 2.5±0.1 (Rv = 3.1), using Baade's Window red 
clump giants. For stellar fields and reddened metal-rich glob- 
ular clusters throughout the Bulge, values of Rv = 3.5 — 3.6 
have been employed (Terndrup 1988, Barbuy et al. 1998). 
We conclude that typical Rv values in Bulge directions can- 
not be significantly lower than Rv = 3.1. 

Assuming that the discrepancy between model and ob- 
served extinction ratio in these strips is entirely due to a 
calibration problem, we estimate a calibration factor of 76% 
to be applied to the Ak,fir values. 

Panels (lib) and (lid) also reveal a strong asymmetry 
between the northern and southern strips. Ak values reach 
up to Ak=1-0 in the north, with most cells in the range 
0.2 < Ak ,2mass < 0.7. In the south, there are few cells with 
Ak,2mass > 0.4. In Sect 4.4, we discuss the north-south 
asymmetry in more detail. 



4.3.2 Strips 3° > |6| > 1° 

As we consider regions closer to the Galactic plane, several 
predictions can be made regarding the extinction values de- 
rived from the 2MASS and DIRBE/IRAS data. First, one 
obviously expects a general increase in both Ak,2mass and 
Ak,fir values. Further, the relation between the two should 
increasingly depart from the identity line, as the contribu- 
tion of background dust is enhanced. Finally, line-of-sight 
variations should also grow in amplitude, as more individual 
dust clouds, with variable dust densities and temperatures, 
are expected to lie along low latitude directions. 

All these predictions are confirmed by inspection of 
Fig. 12, which shows the comparison between Ak ,2mass and 
A k ,fir for the region 3° > b > 1°. The A k ,2mass/A k ,fir 
histograms in panels (12a) and (12c) are dominated by dou- 
ble peak distributions, at Ak,2mass/Ak,fir = 0.75 and 
0.60 for the northern strip, and at Ak,2mass/Ak,fir = 
0.65 and 0.35 for the southern strip. As a consequence, the 
Ak,fir vs. Ak,2mass relation given in panels (12b) and 
(12d) is more complex and not well described by a lin- 
ear fit. Notice that the double peak nature is clearly seen 
in panel (12d) as two distinct branches in the Ak,fir vs. 
Ak,2mass relation. The solid lines shown in these panels 
have angular coefficients derived from the main peak of the 
Ak,2mass/Ak,fir histograms in Panels (12a) and (12c). 

For this region, the model predictions give a foreground 
to total extinction ratio of 81%. If we artificially introduce 
in the model the calibration factor obtained for Ak,fir in 
the previous section, this ratio is reduced to 62%. This is in 
agreement with the typical Ak,2mass/Ak,fir value found 
in the two strips. The main peak in the northern strip (with a 
ratio of 75%) indicates an observed excess of foreground dust 
with respect to the dust exponential model. The secondary 
peak in the north matches quite well the calibrated model 
values. The same applies to the main peak in the southern 
strip, whose ratio is 60%. The secondary peak in the south- 
ern strip (at Ak,2mass/Ak,fir = 35%) suggests that the 
observed Ak,fir is excessively high as compared to the ex- 
pectations of our simple model. This indicates lines of sight 
more strongly affected by dense dust clouds and tempera- 
ture effects caused by dust heated over the DIRBE/IRAS 
temperature correction. 



4.3.3 Strips \b\ < 0.5° 

Very close to the Galactic Plane, the infrared emission by 
hot dust near the Centre, across spiral arms or close to 
star-forming regions, may increase substantially the esti- 
mate of Ak,fir- Therefore, this latter is expected to de- 
viate significantly from Ak,2mass and to display strong 
variations due to scatter in line-of-sight dust. Fig. 13 
shows Ak,fir vs. Ak,2MASS for the |6| < 0.5° strips. The 
Ak,2mass/Ak,fir histograms in panels (13a) and (13c) 
show peaks at Ak,2mass/Ak,fir = 0.27 (27%) for both 
northern and southern strips. The models predict a rel- 
ative extinction of 54% in the same region; introducing 
the Ak.fir, calibration effect found previously, this value 
changes to 41%. Therefore, the observed extinction ratios 
are systematically smaller than the model expectations. As 
discussed in Sect. 4.2, this is just what one expects from the 
existence of non-exponential dust distribution (dense dust 
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clouds) in the foreground and background of the Galac- 
tic Centre, and the overestimation of the Ak.fir, due to 
heated dust. Furthermore, as discussed in the end of Sect. 
2.1, Ak,2mass becomes undetermined for Ak > 2.5 due to 
the limit in the J band optical depth. In fact, this cut-off 
in the Ak,2mass values is clearly seen in panels (13b) and 
(13d), as Ak,2mass saturates for rising Ak,fir values. 



4.4 North-South asymmetry 

A better way of assessing any north-south discrepancy in 
the Ak,2mass and Ak,fir data is to directly compare cells 
with the same value of I but with opposite signs of b. This 
method has the additional advantage that the comparison 
is not affected by the incomplete spatial coverage in the 
southern strip. 

Fig. 14 shows histograms of the ratio of extinction val- 
ues in mirror cells relative to the disc plane. Panels (a), 
(b) and (c) present extinction ratios for mirror cells in the 
strips 3° < |6| < 5°, 1° < |6| < 3°, and 0° < |6| < 
0.5°, respectively. The solid histograms in each panel show 
the distribution of Ak,2m ass {south) / Ak ,2M as s (north) ra- 
tios, while the dashed histograms correspond to the 
AK,FiR(south)/AK,Fin(north) distribution. 

All histograms are normalized by the main peak in 
counts, allowing a better comparison. In the region 3° < 
1 6 1 < 5°, we have two very similar distributions, with peaks 
around ft; 0.6. The general agreement indicates that both ex- 
tinction measurements detect the same asymmetry, with the 
extinction being 60% less in the south than in the north. The 



strong similarity confirms that, in these particular strips, the 
extinction measured in the FIR is dominated by dust on the 
foreground of the stars that were used to measure Ak,2mass- 

In the region 1° < \b\ < 3°, the observed asym- 
metry in the FIR and 2MASS data is again very similar 
in the range < A K , south / A K , north < 0.7, displaying a 
large peak centred around 0.6 (60%). However, the FIR ra- 
tios are systematically larger than the 2MASS one in the 
Ak, south / A K , north > 0.7 domain. This indicates a system- 
atic trend in the sense that southern cells suffer a more sig- 
nificant contribution by background dust than their corre- 
sponding mirror cells in the north. 

In the region |&| < 0.5°, the north-south asymmetry is 
substantially reduced and, interestingly enough, most cases 
now have A K , south / A K , north ~ 1-0. This applies to both 
FIR and 2MASS data, the main difference between the two 
datasets being the larger dispersion in the south/north val- 
ues in the former relative to the later. 

Interpreting this north-south asymmetry is not an easy 
task. The similarity of the histograms in panel (14a) plus 
our model predictions suggest that extinction in the 3° < 
1 6 1 < 5° strips is dominated by foreground structure in the 
form of extended dust clouds (see also map in Fig. 3 and the 
concluding section). At lower latitudes, the extinction values 
increase and become dominated by diffuse dust distribution 
or by dense clouds located close to the Galactic centre. In 
these cases, the shape and position of the Ak, south/ Ak, north 
histograms in both datasets should reflect the global prop- 
erties of the dust distribution rather than the contribution 
of any individual dust cloud. As the FIR data reach deeper 
into the dust columns than the 2MASS, the transition from 
the foreground dominated regime to one dominated by the 
global distribution starts earlier in the former dataset; this 
effect is the likely cause of the difference in the histograms 
in panel (14b). 

In the \b\ < 0.5° strip, the Ak,2mass values are strongly 
limited by optical depth and largely unreliable (Sect. 2.1), 
while the Ak,fir values are suffer from the effects of dust 
heated beyond 21K. Despite these uncertainties the dust 
distribution in this region is reasonably symmetric. There is 
a small trend in the sense of Ak, south/ Ak, north ^ 1? espe- 
cially in the FIR data. This trend is intriguing. The fact that 
most cells have larger Ak,fir values at southern latitudes 
than at northern ones may be telling us something about 
the geometry of the dust distribution relative to the Sun's 
position. 

In order to take into account a possible displacement 
by the Sun relatively to the disk mid-plane, we again use 
our model presented in Sect. 4.1. Considering the Sun dis- 
placed by 5, 15 and 25pc above the Plane, we obtain typical 
south/north ratios in the ts.5 or values of 1.08, 1.27 and 
1.50, respectively. However, the model does not incorporate 
any asymmetry caused by foreground dust clouds, which 
seems to be the dominant effect generating the asymmetry 
in the 2MASS and FIR data. It is thus impossible at the 
present to quantify the Sun's displacement from the disc 
based on the data. We should point out that Unavane et al. 
(1998), using DENIS data, also concluded that the asymme- 
try in the inner Bulge extinction is dominated by foreground 
dust clouds. 
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5 CONCLUDING REMARKS 

We built the Ak extinction map towards the central 10° x 
10° of the Galaxy using the 2MASS Point Source Cata- 
log. We extracted J and K s magnitudes for about 6 mil- 
lion stars in the range 8.0 < K s < 13.0. The adopted map 
resolution is 4' x 4'. It was possible to obtain extinction 
values for ?s80 % of the 32,761 cells defined in the area, 
where 2MASS data were currently available and a Bulge 
giant branch was distinct enough. The extinction affecting 
the bulk of the Bulge stellar population was determined 
by matching the upper giant branch found in the K s , (J- 
K s ) colour magnitude diagram to the reference upper giant 
branch built using de-reddened Bulge fields. The extinction 
values vary from ^4^=0.05 in the edges of the map up to 
Ak=3.2 close to the Galactic centre. The mean extinction 
found is < Ak > = 0.29 with a dispersion a — 0.12; 63 % 
of the cells are within 2-cr of the mean. 

We compared our 2MASS extinction map to that of 
Schultheis et al. (1999) in the region \£\ < 5° and |6| < 1.5°, 
which is common to both studies. Schultheis et al. extinction 
map is based on DENIS photometry. We find an excellent 
agreement between the two extinction determinations, espe- 
cially up to Ak = 1.0. Beyond this limit the values derived 
from the DENIS data are systematically larger. This small 
discrepancy in large extinction regions is not unexpected, 
considering the photometric errors, incompleteness effects, 
and uncertainties in extinction determination. 

We also compared the present extinction map to that 
of Schlegel et al. (1998), which is based on dust emission 
in the far infrared detected by the DIRBE/IRAS instru- 
ments. As the data from the latter are affected by the en- 
tire dust column, with no depth limit, the comparison was 
made separately for regions of decreasing Galactic latitude 
b, which supposedly correspond to increasing contribution 
by dust located on the background of the Galactic Centre. 
The background dust contribution was estimated by means 
of a double exponential dust distribution model. Some ex- 
pected systematic biases, besides that caused by dust on the 
far side of the Galaxy, were also assessed and quantified. 

In general, the extinction values derived from dust emis- 
sion are higher than those from 2MASS, mainly close to the 
Galactic Plane and Centre. We detected two unexpected re- 
gions symmetric and close to the Galactic Centre where the 
two extinction estimates are of the same order. The lack of 
background dust in these low latitude regions could be ex- 
plained by a process of dust grain destruction by UV emis- 
sion from sources associated with continuous star formation 
and /or Post-AGB stars in the central parts of the Galaxy. 

For the cells in the region 3° < |6| < 5°, we observe a 
clear and roughly linear correlation between the Ak values 
from 2MASS and dust emission. We also confirm, as was 
done in Paper I, that the Ak values from 2MASS data are 
in general 73% smaller than those derived from dust emis- 
sion. Since in this region the background dust contribution 
is less than 5 %, the differences between these two quan- 
tities should be smaller than observed. This discrepancy is 
also verified by Arce & Goodman (1999) in the Taurus Dark 
Cloud. It is probably due to systematic effects in the dust 
column density vs. reddening calibration from Schlegel et al. 
(1998), yielding an overestimate of extinction in moderate to 



high extinction regions. We estimate a calibration correction 
factor of 76% for the FIR extinction values. 

For the intermediate 1° < |6| < 3° region, the rela- 
tion between DIRBE/IRAS and 2MASS extinction values 
departs more significantly from the identity line, as expected 
due to the larger contribution by background dust. In this 
region, the typical Ak,2mass /Ak,fir ratio is 65% and could 
be explained by background dust contribution and the cal- 
ibration factor affecting the FIR data. An enhancement in 
the foreground dust with respect to the dust model is ob- 
served in many cells in the northern strip. In the southern 
strip, several cells have Ak,2mass/Ak,fir smaller than ex- 
pected, probably due to dense dust clouds and temperature 
variations, currently not incorporated into our model. 

For the regions very close to the Galactic Plane {\b\ < 
0.5°), we have a typical value for the Ak,2mass / Ak.fir. 
ratio of 27%. Even considering the background dust contri- 
bution and the calibration factor, this ratio is still smaller 
than that predicted by our simple model for the dust dis- 
tribution. This fact is probably due to the overestimation 
of the Ak.fir. values by heated dust above that obtained 
from DIRBE temperature maps. Another possible contri- 
bution to this difference is the existence of systematic ef- 
fects on the Ak,2mass values in high extinction regions 
(Ak,2mass > 2.5), where the 2MASS extinction should be 
significantly underestimated or even unreliable. 

A systematic asymmetry in the Ak values relative to 
the plane of the Galaxy at 1° < \b\ < 5° is observed both in 
the 2MASS and DIRBE/IRAS data. The behaviour and am- 
plitude of this asymmetry with position on the sky suggest 
that the dominant role in creating this north-south asymme- 
try is a more effective presence of foreground dust clouds in 
the northern Galactic strips, such as the Pipe Nebula (Sect. 
2.2). A possible explanation is stellar winds and supernovae 
from nearby OB stellar associations producing dust cloud 
shells (Bhatt 2000). The nearby clouds projected towards 
the central parts of the Galaxy at positive latitudes belong 
to the Ophiuchus dust complex. They are probably related 
to the association ScoOB2, which is at a distance of 145 
pc from the Sun (Bhatt 2000; Onishi et al. 1999). ScoOB2, 
in turn, belongs to Upper Scorpius, which is the eastern- 
most part of the Sco-Cen Association, as studied by means 
of Hipparcos (de Zeeuw et al. 1999) . 

In all regions, significant substructure in the 
Ak,bmassv s. Ak,fir relation is seen, with loops and arms 
stretching out from the main relation. These structures are 
probably caused by intervening dust clouds, with different 
temperatures and densities for different lines of sight. One 
extremely interesting perspective is to model the dust dis- 
tribution within the Galaxy, trying to reproduce as close as 
possible the details of the Ak maps currently available. This 
effort demands models that incorporate, on top of a smooth 
dust distribution, the effects of individual dust clouds, spi- 
ral arms, molecular rings and other structure, possibly with 
variable density contrasts and temperatures. This effort is 
currently under way for the central region of the Galaxy. 
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